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We compare new HERA data for the longitudinal structure function Ef 
with the predictions of different variants of the dipole model. In particular 
we show that the ratio E L /E 2 is well described by the dipole models and 
is rather insensitive to the details of the fit. Fits to E 2 are performed with 
the help of geometrical scaling (GS). Using the property of GS we derive 
the bounds for El/E 2 both for the different versions of the dipole model 
and in the general case. Finally we briefly discuss how the higher Fock 
components of the photon wave function may affect these bounds. 

PACS numbers: 13.85.Ni,12.38.Lg 

1. Introduction 

Recently HI [1] and ZEUS [2] Collaborations have published new data on 
the longitudinal structure function F^(x, Q 2 ) in deep inelastic ep scattering 
(DIS). The HI analysis extends and improves previous data [3], which now 
cover kinematical range from Q 2 = 1.5 GeV 2 /c 2 and x = 0.279 x 10~ 4 up 
to Q 2 = 800 GeV 2 /c 2 and x = 0.0322. ZEUS data has been taken in much 
smaller region from Q 2 = 9 up to Q 2 = 110 GeV 2 /c 2 (see [3j for summary). 
In both data sets there is strong correlation between Q 2 and x values; for 
each Q 2 structure function El (and also F 2 that has been measured in the 
same kinematical points) has been measured over a limited x range, with 
small ads concentrated around small values of Q 2 , see Fig. [lj Moreover, 
since El is difficult to extract experimentally, even recent improved data 
has still large errors. 

Longitudinal structure function is of particular interest since, in the first 
approximation of the parton model, it is equal identically zero [5j (Callan- 
Gross relation) and therefore it is generated entirely by radiative corrections. 


(1) 
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On the contrary, in the dipole model is nonzero, albeit small. Indeed, 
Nachtmann and collaborators have shown that in the dipole model there 
exists a strict bound that El HUH] 

Fl < <7max X F 2 = 0.27 X F 2 . (1) 


This result, hereafter referred to as an EMNS bound, is independent of the 
dipole-proton cross-section, and - strictly speaking - follows solely from the 
properties of the photon-gg wave function. 



Q 2 [GeV 2 ] 


Fig. 1. HI (blue circles) and ZEUS (red 
Fl/Q 2 (right) plotted as functions of Q 
correspond to different Bjorken ads. 



squares) data "1 2] for F 2 /Q 2 (left) and 
|2 . Different points for one value of Q 2 


Using previous HI data [3j the authors of Ref. [8] have shown that the 
bound ([!]) was almost saturated, which is difficult to realize in realistic dipole 
models. In this paper we are going to check if this conclusion remains still 
valid for the new data and what is the value of bound (JT]) for commonly 
used dipole-proton cross-sections. Similar analysis for the Golec-Biernat- 
Wiisthoff model mm has been already performed in Ref. HQ- To simplify 
the analysis we use here the property of geometrical scaling [12] (GS) which 
is to large precision exhibited by the DIS data up to relatively large Bjorken 
ads [13] . We find that for realistic dipole-proton cross-sections the bound is 
indeed lower than 0 with g ~ 0.22 ( i.e 18.5% below the bound of Eq. 0). 
In reality these bounds would be lowered if charm quark mass was included. 

With present experimental accuracy we do not find any significant ten¬ 
sion between Fl/R 2 data and the dipole model. Should such tension arise 
when new data appear, higher order corrections to the dipole model might 
resolve the issue. Therefore we discuss a possibility that corrections to the 
dipole model coming from higher Fock states in the virtual photon may 
change ([!]). Higher Fock components are needed e.g. in the dipole model [9] 
description of the diffractive data flOi . We show that this is possible only if 
there exists a bound for higher Fock components that is analogous to (jT]) , 
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but with longitudinal contribution to F 2 that is significantly different than 
in the case of the qq state. Only explicit calculation of the qqg contribution 
to the virtual photon wave function might give here a definite answer. Such 
calculations have been recently carried out by various authors camaiEi, 
however these results have not been so far applied to the phenomenolog¬ 
ical analysis of DIS. Although the calculation of Ff" / F% 99 with the help 
of Refs. [15, H 6 j might be probably possible, it is beyond the scope of the 
present paper. 

The paper is organized as follows. In Sect.l2] we rederive the EMNS 
bound using geometrical scaling. Next, in Sect. [3pwe fit two versions of the 
dipole model to the present data set for F 2 . To this end we also use the 
property of geometrical scaling. We then compare these fits with the data 
for Tl and discuss fit uncertainties. In Sect. [4] we calculate ratio F\JF 2 
for the afore mentioned fits and compare it with the data and with the 
EMNS bound. We do not find large tension between the data and model 
predictions. An influence of higher Fock states on the EMNS bound is 
discussed in Sect. [5j We conclude in Sect. [ 6 j 

2. Geometrical scaling and the EMNS bound 

For three massless flavors DIS structure functions read HZ]: 

F 2 (x,Q 2 ) = ^3 J dr 2 {|^t(dQ 2 )| 2 + IV’iXdQ 2 )! 2 } a dp (r 2 ), 
Fl(x,Q 2 ) = ^ J dr 2 \iph(r,Q 2 )\ 2 a dp (r 2 ) 
where photon wave functions take the following form 

IVt(d Q 2 )\ 2 = [ dz [^ 2 + (1 - z) 2 ] Q 2 Kf(Qr), 

Jo 

\Mr,Q 2 )\ 2 = 4 / dzz( 1 - z)Q 2 K 2 (Qr) . (2) 

Jo 

Here Ki are modified Bessel functions and 

Q 2 = z{ 1 - z)Q 2 . (3) 

It is convenient to define functions Tt,l 

4>t,l(m = rQ ) = r 2 |^t,l(d Q 2 )| 2 


(4) 
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that depend only on the combined variable u = Qr. Therefore 
F 2 {x,Q 2 ) = ^ j dti{$ T W + ^(a)}^Y^, 

Fl(x,Q 2 ) = ^3 J du T l (r) " dp( ^ /Q) . (5) 


This parametrization is very convenient for the following reasons. First of 
all wave functions •Ft.l^) are universal and do not depend on external 
kinematical variables. Secondly, unlike functions i/’t.lCa Q 2 )> they are ev¬ 
erywhere regular in u. And finally, cross-section (Jdp{u/Q)/u is a localized 
function of variable u that vanishes both for u —> 0 and u —> oo. 

If - as it is in the case of the GBW model - the dipole-proton cross- 
section exhibits geometrical scaling, i.e. Odp(r) = o-dp(^Qs(®)) then the 
integral 

f du d> T L («) adp {U,Q) = [du $T L (u) adp (nQs/Q) = function(r) ( 6 ) 
J u J u 

depends only on a scaling variable 

t = {Q/Qsf ■ ( 7 ) 

Here Q 2 is a saturation scale 

(8> 


Now we can rederive the EMNS bound by considering the ratio 
|r L ( r ) Jdu$ L (u)a dp (u/Q)/u 


F2{T) 


(9) 


J du{$ T (u) + TlCu)} Vd P (u/Q)/u 
and observing that function 


9[U) + 


( 10 ) 


which is plotted in Fig. [2j has a maximum g maiX = 0.2714 for u = 2.591. 
Therefore it follows that 


Fl(x,Q 2 ) 

F2(x,Q 2 ) 


A g'max — 0.27. 


( 11 ) 
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Fig. 2. Solid (red) line: function g(u) defined in Eq. ( |To| . Short dashed (brown) 
line: function g(u) in the case of massive quarks for <5f = mt/Q = 0.5 and long 
dased (blue) line: for S{ = mf/Q = 1.06, which corresponds to the charm quark 
mass m c = 1.3 GeV/c and Q 2 = 1.5 GeV 2 /c 2 . 


The bound (11) has been derived for the case of massless quarks. While 


this is certainly a good approximation for three light flavors, given the fact 
that the lowest photon virtuality in the data set we use is Q 2 = 1.5 GeV 2 /c 2 , 
the inclusion of charm mass effects is going to change ©• For a given flavor 
function g defined in Eq. (10) is no longer a function of scaling variable u only 
but in addition depends on the ratio 5/ = mf/Q 2 . For large Q 2 ( i.e. small 
<5f) g(u,5f) —> g(u). Moreover we have found numerically that everywhere 
in u we have 

g(u,S f ) < g{u) (12) 

and the maxima gmix of g(u,5f) are decreasing with increasing <5f, as illus¬ 
trated in Fig. [ 2 } This is consistent with the observation of Ref. m that 
is a monotonically growing function of Q 2 . Therfore 


0 — 9ma,x(Q ) — Umax ■ 

This allows us to estimate the effect of the charm quark on the ratio: 


(13) 


F, 


(light+c) 


Fl + F, 


(c) 


F l /F 2 + F^ c) /F 2 (c) F 2 (c) /F 2 


(light+c) 


F. 


f 2 + f 2 (c) 1 + f£>/f 2 

1 + gmix/ffmax F^/F 2 


?( C ) 


— 5max 


1 + f 2 (c) /f 2 


ffmax 


(14) 


where the last inequality follows from (13). Note that F 2 . l without any 
superscript refers to the structure function with light flavors only and that 
dependence on Q 2 has been suppressed. We see therefore, in agreement with 
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2-max 

u 0 [l/GeV 2 ] 

A 

x 0 

X 2 /dof 

none 

23.68 

0.389 

0.010497 

1.18 

0.01 

27.11 

0.353 

0.007786 

0.87 

0.005 

29.33 

0.333 

0.006435 

0.79 

0.0005 

38.37 

0.253 

0.003090 

0.70 


Table 1. Parameters of the GBW model fitted to F 2 HI data up to £ max . 


Ref. [8], that inclusion of charm lowers the bound on ' s c /F 2 ) " 1 ' m 
proportion that depends on /F^- For m c = 1.3 GeV/c and for the lowest 
Q 2 in the present data set we get numerically g j c L ~ 0.05 (see Fig. Q), 
which gives giiL/gmax ~ 0.19. We do not know what is the fraction of charm 
in the present data sample, however taking a typical value of _F 2 /F? ~ 25%, 

we get that gS t+c ' ~ 0.23. This means that bound (11) is lowered for the 
lowest Q 2 bin by approximately 18% and is approaching for higher Q 2 . 
In the following we will ignore charm contribution and stick to the bound 
(fTTj) coming back to this point in Sect. [4} 


3. Dipole models and geometrical scaling 

In order to check how far the bound ( p~T| ) overestimates actual predictions 
of the dipole model with realistic dipole-proton cross-section, we are going 
to compute ratio 0 for a given erdp in terms of scaling variable r. We 
will see that for commonly used parametrizations of <7dp> ratio F\J Fz is a 
slowly varying function of t with a maximum equal approximately 0.216 — 
0.224, which only slightly depends on the parametrization actually used. 
To this end we have decided to perform our own fits to the F 2 data over 
the restricted kinematical range where the longitudinal structure function 
Fj J has been measured by HI. The reason for this is threefold. Firstly, new 
data is of much better quality than the previous data to which the dipole 
model parameters have been fitted. Secondly, we do not aim at a global fit, 
but rather at a fit which covers only the points where Fp has been measured 
as well. Therefore fit parameters - as we shall see in the following - will 
be different from the ones obtained in the global fits. And finally, we have 
decided to fit the data with the help of geometrical scaling - a procedure 
not used so far in the fits to the DIS data. 

Fitting dipole models to the data becomes very easy when F 2 depends 
only upon single scaling variable r. This happens because points corre¬ 
sponding to one particular value of Q 2 but different x’s (see ’’stacks” in 
Fig. 0 correspond to different values of r and are therefore shifted horizon¬ 
tally - if plotted in terms of t - by values that are different for different 
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x’s. As a consequence dipole model predictions fall on a universal curve 
(up to an overall normalization a o), and data fitting consists in changing do 
and the parameters defining scaling variable r, i.e. xq and A. By varying 
these three parameters one forces experimental points to fall on theoretical 
prediction, rather than by changing theoretical predictions one is trying to 
reproduce experimental points. Therefore this method is very efficient, as 
it does not require time consuming recalculations of the theoretical curve. 

Even in the case of dipole-proton cross-sections that violate GS by ex¬ 
plicit (albeit weak) dependence on x, like in the model of Iancu, Itakura 
and Munier P33, it is still possible to make a GS-like fit by keeping x at 
some fixed average value x ave and then study the uncertainty of theoretical 
predictions by changing x over the range that is covered by experimental 
data. We shall come back to this point later. 

Let us first consider the simplest version of the dipole model, namely 
the GBW parametrization [9j, which - up to an overall normalization do - 
takes the following form in terms of scaling variable t 


r GBW / 
J d P 


FXP"{u,t) ^ f u 2 \ 






(15) 


Pluging (15) into Eqs.([5]) gives unnormalized theoretical prediction for the 
structure functions divided by Q 2 , which will be denoted in the following by 
small characters fi Experimental data are tabulated in a set of discrete 
points {Qi,Xi}, and we fit three free parameters of the model, do, xo and A, 
by transforming experimental entries in the following way: 


F 2 {xi,Q 2 ) 


Q 2 a o 


F 2 (x i ,Q 2 ) = f 2 


Ti = 


Q 2 


Qi(xi 


(16) 


and demanding that they are equal to the theoretical prediction at the 
pertinent value of scaling variable r, with Q s given by Eq.Q. The results 
are shown in Table [lj Since GS is supposed to work the best for small values 
of Bjorken x’s we have performed a number of fits restricting the allowed x 
region up to a maximal value denoted as x max . We see that even without 
any cut on the maximal value of x, i.e. for x as large as 0.0322 (the highest 
x in the analyzed data set) \ 2 of the fit is quite reasonable. By restricting 
analyzed data to the smaller and smaller range of Bjorken x’s y 2 is getting 
smaller, but also model parameters vary substantially. Parameter do is 
much smaller than in the original fit of Ref. m d£ BW = 23 mb = 59 GeV" 2 . 
Exponent A approaches the value of Ref. [9] A GBVV = 0.288 only for small 
x max (note that maximal x in Ref. |9] was equal to 0.01, whereas the lowest 
x = 6 X 10~ 6 was much below the minimal x of present analysis), otherwise 
being consistent with model-independent analysis of Ref. m- The results 
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of the fits, together with the original parametrization of Ref. |9J are plotted 
in Fig. [3j One can see rather good agreement of fits from Table [l] with the 
data, and - also quite importantly - good quality of GS of the data. 



Fig. 3. HI (blue circles) and ZEUS (red squares) data [TJ [2J for F 2 /Q 2 plotted as 
a function of scaling variable r for different values of fit parameters corresponding 
to the GBW model given in Table [l] Fit parameters in the first panel correspond 
to the original fit of Ref. |9j with no charm data included. 


Finally in Fig. [4] we plot data for Fl/Q 2 as a function of t together with 
theoretical parametrization of Ref. 0J and the present fits corresponding 
to Tab. [lj We can see that due to still large experimental errors of F\ J all 
parametrizations, although different, describe well the data. 

As a second example let us consider a dipole model by Iancu, Itakura 
and Munier (IIM) |18j . where the dipole-proton cross-section is defined in 
terms of two functions 


u 


ln(4r/« ) 

-\ !' 2k\\h(1/x) 


Ai(u,t,x) = A 0 J 
A 2 (u,t) = 1 — exp (— a\n 2 (bu/y/r)) 


( 17 ) 
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Fig. 4. HI (blue circles) and ZEUS (red squares) data [HE: for iff/Q 2 plotted as 
a function of scaling variable r. Solid curves correspond to the GBW model fits of 
Table [lj 


and 



A\(u,t, x) for u 2 < 4t 
A2(u,t) for 4 t<u 2 


(18) 


Here 7 = 0.63 is related to the anomalous dimension of the forward scat¬ 
tering amplitude in the BFKL formalism, while n = 9.9 corresponds to the 
diffusion coefficient. Parameters a and b are determined uniquely by gluing 
Ai and A 2 and their derivatives at u 2 = At. Parameter Aq is in princi¬ 
ple free, however, as it was shown in Ref. [18] the best \ 2 was obtained for 
Aq = 0.7 and for the purpose of the present work we will keep it fixed at this 
value. Therefore the only free parameters are, as in the case of the GBW 
model, an overall normalization < 7 o and two parameters of the saturation 
scale: A and xo- 

However, there are two important differences between IIM and GBW 
parametrizations. First of all, for small value of u amplitude A\ exhibits 
explicit violation of GS, since it depends both on r and x. For the purpose 
of the present work we will keep x entering A± fixed at the average value 
x ave calculated for the interval where the fit is performed. The accuracy of 
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^max 

^ave 

u 0 [l/GeV 2 ] 

A 

x 0 

y 2 /d°f 

none 

0.00359 

20.22 

0.597 

0.002553 

1.76 

0.01 

0.00182 

21.50 

0.583 

0.002140 

1.57 

0.005 

0.00121 

25.56 

0.531 

0.001392 

1.31 

0.0005 

0.00022 

34.30 

0.389 

0.000645 

0.75 


Table 2. Parameters of the IIM [TS] model fitted to F 2 HI data up to a; max . 






Fig. 5. HI (blue circles) and ZEUS (red squares) data [H'^i for F 2 /Q 2 plotted as a 
function of scaling variable r for different values of IIM model fit parameters given 
in Table [2] 


this procedure is then checked by putting in (18) x equal to the maximal 
and minimal value of x in a given interval. It will turn out that the structure 
functions are sensitive to this variation of x at the level of a few percent, 
however the ratio F\J F 2 is almost independent. Next difference appears 
due to the dependence of A\ on A. To solve this issue we have adopted 
an iterative procedure, consisting in fixing A at some initial value, then 
performing the fit and plugging in the fitted value of A to Ai and performing 
the fit again. Usually after four, five steps satisfactory convergence has been 
achieved. The results are given in Table [2} 

By inspecting Table [2] we see that the quality of fits is worse than in the 
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Fig. 6. F 2 (t)/Q 2 for the parameters from the third row of Table [ 2 ] with x in A\ 
from Eq. ( p~7| ) replaced by :r m in (upper curve), x ave (solid middle curve) and x max 
(lower curve). 


case of the GBW model. This is in contrast with the original fits of Ref. ||18] 
which, however, were performed over the data set covering much lower x’s 
than in our present analysis. One should also note that the errors of the 
old data sets are bigger than the ones of the present data. This is also the 
reason why fit parameters are different than in our case. For illustration 
purposes we have plotted in Fig.JbJi^M/Q 2 for the original IIM parameters 
and for three choices of x max from Table [2} Magnifying first plot in Fig. [5] 
one could see that for r > 10 the original IIM curve missed the experimental 
points, which have rather small errors. 

In order to check sensitivity of the IIM fits to the fact that x dependen- 
dent piece of A\ amplitude © has been replaced by a constant value ^Cave? 
we plot in Fig. [fi] E 2 Cr)/Q 2 f° r the parameters from the third row of Table [ 2 ] 
with x in A\ replaced by x m i n (upper curve) and x max (lower curve). For 
better resolution the plot is restricted to r < 5. We can see that theoretical 
uncertainty introduced by this procedure is in fact much smaller than the 
experimental errors. 

Finally in Fig. [7] we plot F^/Q 2 as a function of r for three sets of 
parameters from Tableland for the original set of parameters from Ref. |18| . 
One can see that all curves describe the data reasonably well due to the large 
error bars of F\,. 


4 . EMNS bound for dipole models 

Now we are able to compare the EMNS bound with F\JF 2 ratio calcu¬ 
lated in the dipole model for realistic dipole-proton cross-sections <7d p dis¬ 
cussed in the previous Section. One should note that for the GBW dipole 
model F l /F 2 does not depend on the values of xq, A and erg. It is not the 
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Fig. 7. HI (blue circles) and ZEUS (red squares) data [TJ [2] for F^/Q 2 plotted as 
a function of scaling variable r. Solid curves correspond to the IIM model fits of 

Table [2 



Fig. 8. Ratios FjJF± as functions of r for the GBW dipole model and for the IIM 
model with parameters corresponding to the fit with x < 0.005. The scale of the 
plot has been enlarged in order to make small differences between the two curves 
visible. 


case for the IIM dipole model, but we have checked explicitly that for all 
parametrizations of Table [2] the differences in Fj J / F 2 are negligible. There¬ 
fore in Fig. [8] we plot ratios F\JF± as functions of r for the GBW dipole 
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model 

7m ax 

(Fl/F^ m&x 

GBW any 

1.165 

0.224 

IIM all x 

1.417 

0.217 

IIM x < 0.01 

1.411 

0.216 

IIM x < 0.005 

1.413 

0.216 

IIM x < 0.0005 

1.418 

0.217 


Table 3. Maxima of FjJ F 2 for different fits to F 2 . 


model and for the IIM model with parameters corresponding to the fit with 
x < 0.005. We see that ratios F\JF 2 are in fact almost model-independent. 
This is further confirmed in Table [3] where we collect the maximal value of 
Fj J / F ‘2 for the GBW and IIM parametrizations. 

One should note that each curve in Fig. [8] corresponds to a different 
definition of scaling variable r, so one cannot superimpose experimental 
data on that plot. This is done in Fig. [9] where we plot F\J F 2 for the 
unrestricted fits of the GBW and IIM models corresponding to the first rows 
of Table [l] and [2] respectively. The errors of the ratio have been calculated 
neglecting correlation between errors of F 2 and Fl- 


A 



Fl 

f 2 


/ 


A F l 
F l 


2 


+ 



(19) 


This procedure overestimates the errors, however, given the fact that Fl < 
0.27F 2 and that experimentally absolute errors of Fl are 2-10 times larger 
than AF 2 , the error of the ratio (19) is determined to very high precision 
by AFl alone. We shall come back to this point later. 



T T 


Fig. 9. Ratios Fl/F 2 plotted as functions of scaling variable r . Straight line 
corresponds to the EMNS bound. Solid line in the left panel corresponds to the 
unrestricted GBW fit of Tab. [I] (first row), and in the right panel to the unrestricted 
IIM fit of Tab. [2] (first row). 
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We can see from Figs. [8] and [9] that for realistic <7d p theoretical predictions 
lie below the EMNS bound. Indeed, we see that the maximum of F\ J /F 2 
is of the order 0.216 — 0.224 and only slightly varies from fit to fit. This 
is illustrated in Table [3} A question arises whether data points - which 
for some values of r exceed the EMNS bound - are indeed, as suggested 
by the authors of Ref. [8], saturating bound s being as a consequence 
incompatible with the dipole model. To this end we have simply calculated 
X 2 of F h /F 2 for the unrestricted GBW fit and obtained very good result: 
0.7425. To check whether this value is affected by the fact that we have 
not taken into account correlations between errors of El and F 2 , we have 
recalculated x 2 neglecting F 2 errors, which gives x 2 that changes by less 
than 1%. Therefore indeed, as already mentioned above, x 2 value is totally 
driven by the large errors of El. 

It is interesting to check whether an overall shift of the dipole model 
prediction for F\ J /F 2 would improve agreement with the data. In this way 
we shall have a quantitative measure of the quality of the dipole model 
prediction for the ratio F\JF 2 and also an indication how much room is 
there for the higher order corrections that we are going to discuss in Sect. [5] 
To this end we use the GBW fits allowing for arbitrary normalization of the 
ratio 


El 

F 2 



( 20 ) 


and calculate x 2 /d.o.f. (assuming one degree of freedom, namely Af) as a 
function of Af. This is illustrated in Fig. [TO] for two GBW fits of Table [l] 
unrestricted x and x < 0.0005. We see that depending on the fit the data 
prefers Af sligthly smaller (unrestricted x) or sligthly larger [x < 0.0005) 
than 1. These are negligible changes and therefore at this stage we conclude 
that the data for the ratio El/ F 2 is compatible with the dipole model. One 
should perhaps remind here again, that this ratio is quite stable, even for 
fits to F- 2 that are visibly different. 

One can also see that lowering Af by 20%, which would be requried by 
taking into account charm mass effects ( c.f .. discussion at the end of Sect. [2]) 
increases x 2 , but only by a small ammount. Of course the detialed study 
of charm mass effects would require to go beyond GS used in this analysis 
- which in turn would result not only in the change of Af, but also in the 
change of the shape of F\ J /F 2 - but no drastic difference to the present 
analysis should be expected. 


5. EMNS bound and higher Fock states 

In the dipole model the virtual photon dissociates into a qq pair which 
subsequently interacts with the proton target. However, it is clear that 
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Fig. 10. Change of \ 2 of A f x Fl/F-j as a function of TV" for the GBW fits from 
Table [I] for unrestricted x (left) and x < 0.0005 (right). 


higher Fock states have to contribute as well, similarly to the diffractive 
DIS where the next Fock component, namely the qqg state, is dominant at 
small j3 [10]. Full calculations of the qqg component of the photon wave 
function have recently appeared in the literature nanansi, however they 
have not been yet applied to the phenomenological analysis of DIS. 
Structure functions in the dipole model are given as an expansion 

F 2tL (x,Q 2 ) = F^\x,Q 2 ) + F^f\x,Q 2 ) + ... . (21) 

The EMNS bound derived in Sect. [2] is in fact valid only for the first 
component of © 


G 


Fl qq \x,Q 2 ) 
F!f q \x,Q 2 ) 


^ ^max 


0.27. 


( 22 ) 


One should note, however, that at this order of perturbative expansion the 
loop corrections to the leading order Fock component may change the value 
of the bound (22). Similarly for the next Fock component we would have: 


Fj qq9 \x,Q 2 ) 
F^ qq9) (x,Q 2 ) 


< hr. 


(23) 


where 


0 f; Vax T 1- 


(24) 
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Up to this order one can derive the modified EMNS bound: 
Fl(x,Q 2 ) 


F^ qq \x,Q 2 ) + F^ qq9 \x,Q 2 ) 


F 2 (x, Q 2 ) F.^ q \x, Q 2 ) + F^ qq9 \x, Q 2 ) 

Smax F^ qq) (x, Q 2 ) + h m ax F(x, Q 2 ) 


< 


F 


(<?<?)/ 


2 (uQ 2 ) + i ? 2 


(mg ), 


[x, Q 2 


— 9 max .' 


1 + 5s(x, Q 2 ) 
1 + s(x, Q 2 ) 


(25) 


where 


and 


<5 = 


0 < 5 < 


~ 3.7 


5max 


e(x,Q 2 ) = 


9 max 


F^ qq9 \x,Q 2 ) 

F?\x,Q 2 ) 


(26) 


Since the dipole model with qq component only describes F 2 rather well 
over the wide kinematical range, we do not expect e to be large. For the 
purpose of the present analysis we assume that e does not exceed 20 %. 



Fig. 11. Modification factor / of Eq. (27) as a function of e and S. Light green 
plane corresponds to / = 1. 

In order to change the value of the EMNS bound we need the modifica¬ 
tion factor 


1 + £ 


(27) 
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to be significantly different from 1. In Fig. 11 we plot f(e, <5) for 0 < e < 0.2 
and 0 < 5 < 1 /g max . We see from Fig. 11 and from Eq. (|27j) that / = 1 for 


<5=1 and that / gets smaller than 1 if 5 < 1 and / > 1 for 5 > 1, and that 
the difference | f — 1| is growing with e. From Fig. 11 one may conclude that 
the modification of the EMNS bound by more than ±20% might be rather 
difficult. Most probably higher Fock components would modify by less 
than 10%, but to quantify this statement one needs to calculate explicitly 
limax which is beyond the scope of the present paper. One should note at 
this point that in this case an inclusion of mass corrections due to the charm 
and possibly bottom quarks should be included, as these corrections would 
be of the same order or even larger than a contribution from the higher Fock 
components. 

To conclude this Section let us only remark that e depends on kinemat- 
ical variables, so one cannot exclude a’priori a situation that there exists a 
kinematical corner where a correction to (11) is of importance. 


6. Conclusions 

Dipole model offers effective and intuitive description of deep inelastic 
scattering which goes beyond leading twist approximation. In this paper 
we have used different variants of the dipole model, i.e. two different forms 
of photon-proton cross-section fitted to ±2 data over different ranges of 
Bjorken x’s. We have decided to perform our new fits in the kinematical 
range where both F-i and Fl have been measured. This kinematical region 
does not extend to very small x’s as it is in the case of the recently published 
combined HERA data m and therefore fitted parameters are different than 
the ones of global fits. For the same reason the IIM model [18] that has 
been specifically devised for low x region gives larger y 2 than the simplest 
version of the the GBW model j9j, which quite satisfactorily describes F?. 

In order to fit model parameters we have used the property of geomet¬ 
rical scaling, which boils down to the fact that data points of the same Q 2 
but different x’s disperse when plotted in terms of scaling variable r and 
fall on one line (compare Figs. [I] and e.g. [3]). We have found that GS is 
present both in the case of F% and Fl as well. One should, however, take 
this property with care for the present set of data, since for each value of Q 2 
only a few points of different Bjorken x’s have been measured. This x — Q 2 
correlation of the DIS data, particularly pronounced in the present case of 
Fl, is a common obstacle in deriving firm conclusions on the quality of GS. 
In this context it is also worth mentioning that there exist different sources 
of GS. The first one is related to the genuine property of the initial state 
of the proton target described theoretically by the Color Glass Condensate 
formalism [20] and nonlinearities in the parton evolution m, which are of 
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importance at small values of Bjorken x. The other one is a property of the 
linear parton evolution equations [22] , which build GS at larger Q 2 and at 
not necessarily very small x’s. 

We have next compared dipole model predictions - with parameters 
fixed by fits to F 2 - with F l data. Here, due to large experimental er¬ 
rors, agreement is quite good. Next we have studied ratio of the structure 
functions G = F\JF 2 for which in the dipole model there exists a strict 
bound y = 0.27 0 derived in Ref. [8j. We have rederived 0 with the 
help of geometrical scaling. We have also shown that in the dipole models 
discussed above gmax ~ 0.216 — 0.224 which is approximately 18.5% below 
the EMNS bound. Different fits give very similar ratio g , which is only 
residually dependent on the values of fit parameters (this concerns only the 
IIM dipole model, since ratio F^/Fz is parameter independent in the GBW 
case). Comparing G = Tl/T 2 with the data we have established that the 
GBW model reproduces G with high precision. We do not see any tension 
between the data and the dipole model as far as ratio G is concerned, even if 
charm mass effects, which lower the EMNS bound, are taken into account. 

Dipole model and the EMNS bound discussed so far rely on the first Fock 
component of the photon wave function. Including higher Fock components, 
like a qqg state, might in principle change theoretical prediction for q mM . 
For this to happen, longitudinal part of the qqg state compared to the 
transverse one has to be significantly different than in the qq case. 
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